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We present a new theoretical framework for describing an impurity in a trapped Bose system in 
one spatial dimension. The theory handles any external confinement, arbitrary mass ratios, and a 
weak interaction may be inclnded between the Bose particles. To demonstrate our technique, we 
calculate the ground state energy and properties of a sample system with eight bosons and find an 
excellent agreement with numerically exact results. Our theory can thus provide definite predictions 
for experiments in cold atomic gases. 


PACS numbers: 03.75.Mn,67.85.-d,68.65.-k,67.85.Pq 


and a single scat- 
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An impurity interacting with a reservoir of quantum par¬ 
ticles is an essential problem of fundamental physics. Fa¬ 
mous examples include a single charge in a polarizable 
environment, the Landau-Pekar polaron Qi , a neutral 
particle in superfluid "‘He Q, a magnetic impurity in a 
metal resulting in the Kondo effect 
tering potential inside an ideal Fermi gas [^, 
latter system is famous for the Anderson’s orthogonality 
catastrophe Q- In these settings the impurity behavior 
can provide key insights into the many-body physics and 
guide our understanding of more general setups. 

A complicating feature of many impurity problems is 
the presence of interactions at a level that often pre¬ 
cludes the use of perturbative analysis and self-consistent 
mean-field approximations. This implies that analyti¬ 
cal approaches are highly desirable and exact solutions 
are, when available, coveted tools for benchmarking other 
techniques. This is particularly true for one-dimensional 
(ID) homogeneous systems where solutions can often be 
found based on the Bethe ansatz [8-13|. These solu¬ 
tions are the essential ingredients for our analytical un- 
derstandingof highly controllable experiments with cold 
atoms For instance, the exactly solvable prob¬ 

lem of the single impurity in a ID Fermi sea 0 - the 
Fermi polaron - can be used to study the atom-by-atom 
formation of a ID Fermi sea (^ . 

While Fermi polarons have been studied intensively in 
recent times using cold atomic setups both experimen¬ 
tally and theoretically HI-2^, the physics of impurities 
in a bosonic environment is only now becoming a frontier 
in cold atom experiments |26l - l29l| . This pursuit requires 
theoretical models for describing the Bose polaron [ 13 - 
Idsj l. where, in contrast to the Fermi polaron, an exact 
solution is not known even for a homogeneous ID sys¬ 
tem. Here we provide a new theoretical framework that 
captures the properties of an impurity in a bosonic bath 
confined in one spatial dimension. Our (semi)-analytical 
theory thus provides a state-of-the-art tool for exploring 
the properties of Bose polarons in ID. 

The proposed framework works with a zero-range po¬ 
tential of any strength and handles any number of ma¬ 





FIG. 1: (color online) a) A sketch of the doubly-degenerate 
ground state for the system of one A and eight B parti¬ 
cles trapped in a one-dimensional harmonic potential for in¬ 
finitely strong repulsive inter-species interaction. The ± cor¬ 
responds to the parity of the system. Panels b)-d) show the 
AB pair-correlation function: the expectation value of the 
S(za — x)5{zb — y) operator, where Za{b) is the coordinate of 
particle A{B). The system under consideration is the same as 
in a) although with finite inter-species interaction strength, g. 
Panels e)-g) show the pair-correlation function for different 
intra-species interaction strength gsB for fixed g. 


jority particles in external confinement of various shapes 
which is beyond known analytical approaches to this 
problem. Our method is also applicable to describe ex¬ 
perimental setups that have different trapping potentials 
for the impurity and majority particles and different mass 
ratios. In addition, weak majority interactions may be 
included using the well-known Gross-Pitaevskii equation 
(GPE). While our method is not exact, we have bench- 
marked the energetics and density profiles against numer¬ 
ical results (d^ l and find agreement for up to ten particles 
at the level of a few percent. To illustrate the method in 
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this letter, we examine a system of eight bosons and an 
impurity in a harmonic trap. Fig. HJi) shows a sketch of 
this system with vanishing boson-boson and large boson- 
impurity interaction. This leads to separation of the two 
components. Notice that the ground state must retain 
parity and is thus a linear superposition of the two spa¬ 
tial configurations outlined. Using the pair-correlation 
function we can clearly see in Fig. [D3)-d) how the im¬ 
purity moves to the edge of the system as a function of 
the inter-species interaction. Increasing the intra-species 
interaction we witness the opposite effect as the impurity 
goes to the center of the system as seen in Fig. [I^)-g). 

Formalism Our two component system consists of 
one type A (impurity), and Nb identical type B bosons 
(majority) with masses mA and m b respectively. For the 
sake of argument, we confine particles in harmonic poten¬ 
tials with trapping frequency ujb for the bosons and uja 
for the impurity. In this letter we adopt harmonic os¬ 
cillator units for the majority particles, i.e. we measure 
length in units of 6 = yJ’hjmB'^B and energy in units 
of /kus. Accordingly, the Hamiltonians, for the impu¬ 
rity atom with coordinate x and a majority atom with 
coordinate ?/, are expressed as 


Ha{x) 


pI + m^B^ABX^ 

2mAB 


, Hb^v) 



( 1 ) 


where mAB = ujab = and p denotes 

the corresponding momenta. The interaction between A 
and B particles is assumed to be of a short range and 
hence modeled by the Dirac delta-function with strength 
g. The boson-boson interaction is also given in the stan¬ 
dard pseudopotential interaction model with cou¬ 
pling constant qbb ■ Both interaction strengths are given 
in units of [bhujB]- The overall Hamiltonian of the sys¬ 
tem is iJ = Ha{x) -f HBivi) + gS{x-y^) + 

J2i<k PBB^iVi — Vk), where yi are the coordinates of the 
bosons (see Supplemental Material 0). 

In order to find the eigenspectrum of the Hamiltonian 
for arbitrary g we introduce a new (semi)-analytical ap¬ 
proach. More specifically, we consider the impurity as 
the ‘slow’ variable and introduce the adiabatic decompo¬ 
sition of the total wave function 


^{x,yi,...,yNB) = ^4>o{x)^3i.yi^----:yNB\x), ( 2 ) 

i=i 

where is the jth normalized eigenstate of the eigen¬ 
value problem HBiyi)^j = Ej{x)^j which we solve 
assuming that the impurity gives rise to a zero-range po¬ 
tential at a fixed position x (see Supplemental Material 
0). First we consider an ideal Bose gas, i.e. gBB = 0, 
and write (yi,..., j/atb |x) = f^{yi\x) with a 

symmetrization operator S (acting on the yi coordinates) 
and fi. 3 {yi\x) being the kfth normalized eigenstate of 
HB^yi) for a given x. Notice that every function {yi\x) 



FIG. 2: (color online) The energy spectrum of low-lying states 
for a system of eight non-interacting bosons and an impurity 
with mAB = kJAB = 1. The inset shows the derivative of 
the ground state energy dE/d{ — l/g) divided by Nb, which 
is related to Tan’s contact parameter. The lines are guides 
for the eye. 


has a discontinuous derivative at yi = x due to the zero- 
range interaction. This is quantified with the standard 
delta function boundary condition which dictates that 
the difference in the slopes of the wave function, from 
the left and right sides of x, times a \/{2g) factor must 
be equal to the value of the wave function taken at yi = x. 
As 1/g —>■ 0 the wave function must therefore vanish at 


yi = X. 

We can include interactions among the majority par¬ 
ticles under the assumption that these may be described 
by the ID GPE (see Supplemental Material 0 |). It has 
been previously discussed that this is an accurate de¬ 
scription for weak interactions between the bosons 51- 
0 . In this case we need to use a dressed single-particle 
wave function /^j instead of fj^j{yi\x). The function fy 
satisfies the ID GPE complemented with the boundary 
condition at yi = x, 


P{x)fy = + ^yf + ■ OBBlfki P ) /fc^ (3) 


where yt{x) is a chemical potential, gBB is determined 
through the three-dimensional boson-boson scattering 
length, Os, as gBB = ^ ’ where ljbi and 

are the two frequencies in the directions of strong confine¬ 
ment 5^ . The boson-impurity coupling constant can be 
also related to the corresponding scattering length |55|. 
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After determining the function we obtain a coupled 
system of equations for (j)j (x) (see Supplemental Material 
[49l|l. The coupling terms in this system correspond to 
the transition of a majority particle from fk to fk' ■ For 
bosons, this is a coherent process which contributes sig¬ 
nificantly if the impurity is placed in a region with high 
density of majority particles. Physical intuition, how¬ 
ever, tells us that in the ground and low-lying excited 
states the impurity is pushed to the edge of the trap if 
gsB ^ 9 and Nb ^ 1- Otherwise, the impurity would 
deplete the majority particles notably from the ground 
state of the one-body harmonic oscillator which is very 
expensive energy wise for Nb ^ Hence, for large Nb 
we neglect the coupling terms between different (f>j{x), 
which rigorously gives us an upper bound for the ex¬ 
act energy of the ground state [^. However, we expect 
this approximation to be very accurate also for low-lying 
excited states. We can therefore obtain (j)j{x) by solv¬ 
ing numerically a single differential equation. From the 
mathematical point of view the presented approach is 
similar to the Born-Oppenheimer approximation or the 
hyperspherical adiabatic method [57|. The physics is, 
however, different. Indeed, we develop our method for 
a many-body system where we expect that for the same 
computational time the relative precision is increasing 
with the number of particles. Clearly, the discussion 
above applies to a Bose polaron problem in arbitrary 
confinement. Moreover the trapping potentials as well 
as the masses can be different for the A and the B parti¬ 
cles. Below we show that our framework compares very 
well to exact numerical results for the equal mass case 
mA = niB- In the heavy impurity case with niA > tub, 
we expect our model to perform equally well as the im¬ 
purity becomes increasingly stationary. The final case of 
a light impurity niA < btib will not be studied further 
here, but we note that since it must also go to the edge 
of the trap, the arguments above for neglegting couplings 
terms still hold, and we expect our model to work well 
also in this case. 

To conclude the presentation of our method we com¬ 
pare its predictions with the exact results obtained us¬ 
ing the numerical approach developed in Refs, [i^, 
for 9 bb = OjiBAB = ^AB = 1- We find that the rel¬ 
ative precision of the method increases with Nb and 
we pick a sample system with Nb = 8. We start by 
analyzing the energies in Fig. [2] Our model yields re¬ 
sults that are slightly above the numerically exact values 
with a maximum deviation of a few percent. Next we 
check that the model reproduces the derivative of the 
ground state energy with respect to the coupling con¬ 
stant, dE/d{—l/ 9 )/NB, vs. Nb, see the inset in Fig. [2] 
This derivative for fixed g determines the probab ility 
for a given boson to be close to the impurity [H, 1^ . 
We see that for a large number of bosons this proba¬ 
bility becomes smaller, manifesting that the impurity is 










x/b y/h 


FIG. 3: (color online) Panels a)-b) show the density dis¬ 
tributions for the impurity, a), and majority, b), atoms for 
Nb = 8 and different values of the inter-species interaction 
strength, g. The numerically exact results are shown with 
blue dots in a) for the three corresponding values of g, and 
in b) for g = 50. Notice that the label in a) also applies to 
b) and vice versa. c)-d) The density distributions for a fixed 
value of p = 1 but different intra-species interaction strength, 
gsB- e)-f) The density distributions for different values of 
mass ratios, rriAB ~ 1, 3, and 10 with gBB = 0. 


pushed far from the center of the trap. It is also inter¬ 
esting to note that for large Nb we find numerically that 
dE/d{—l/g) is almost independent of Nb- 

We have also compared density profiles and pair- 
correlation functions for the ground state and again find 
only minute differences, see Fig. and b). Note that 
the impurity density splits for large g. On the other 
hand, the majority particles are almost unperturbed by 
the interaction. This means that an adiabatically slow in¬ 
crease of g moves the impurity to the edge of the system, 
as also shown in Fig. [I}i)-d). As the number of bosons 
increases the impurity gets pushed further towards the 
edge of the trap, and dE/d{—l/g) /Nb decreases. In the 
energy domain it leads to a doubly degenerate ground 
state at 1/g = Q since the impurity can be pushed to ei¬ 
ther the left or the right edges of the trap, see Fig. [2]and 
Fig. [Di). This should be contrasted with the Fermi po¬ 
laron system where for large interaction the ground state 
is Nb + 1 times degenerate and the impurity is localized 
in the middle of the trap [H^, El, 621. 


Results To further illustrate the model we continue 
our discussion of an impurity interacting with eight 
bosons. However now we allow weak interaction between 
the bosons as well as different mass (or frequency) ratios. 
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tion is almost a Gaussian function, and therefore the mo¬ 
mentum distribution will also assume a Gaussian form. 
Notice the characteristic oscillations in the wings of the 
distributions which could be very helpful for the exper¬ 
imental detection of the Bose polaron. The momentum 
distribution for majority particles is not plotted because 
there is no noticeable change in the distribution as we 
change g and/or itiab- 


FIG. 4: Snapshots of the momentum distribution for the im¬ 
purity with different values of g for itiab = 1 (left) and dif¬ 
ferent values of btab for g — 1 (right). We assume pss = 0 
and UJAB = 1. 


For large interactions this setup is already beyond cur¬ 
rent numerical approaches. We stress once again that in 
our method the numerical complexity does not increase 
with Nb^ and Nb = 8 is chosen as a particular example. 


First we fix g = l,mAB = = 1 and change 

the intra-species interaction strength, gsB- In Fig. [3}:) 
and d) we show the density for the impurity and major¬ 
ity atoms in the ground state. In this case the energy 
is minimized if the impurity is pushed towards the mid¬ 
dle of the trap. This is readily understood for the case 
g = 9 bb where the impurity particle should have the 
same density distribution due to the boson-impurity ex¬ 
change symmetry of the Hamiltonian. In our case the 
difference in the densities of majority and impurity for 
g = gBB = 1 is due to the different treatments of these 
components. Notice also that we expect the GPE and our 
approach to yield only qualitative results for such a large 
boson-boson interaction strength. It is also worth notic¬ 
ing that if both interactions g and gBB are very large one 
can approach the problem directly using the Bose-Fermi 
mapping , where one also expects enhancement of the 
impurity density in the middle of the trap compared to 
the = 0 case. Next we compute the pair-correlation 
function which again demonstrates that the impurity is 
situated close to the origin, see Fig.[T^)-g). Consider now 
different masses for A and B particles for = 0,g = 1 
and UJAB = 1 [S^- As shown in Fig. [3]e)-f), when btab 
becomes larger the external potential localizes the impu¬ 
rity in the middle of the system. For ttiab oo the 
impurity constitutes a delta-function barrier in the mid¬ 
dle of the harmonic trap, the solution to which can be 
found in Ref. 63 • From this picture it is apparent that 
the density of the majority particles should be suppressed 
at the origin, see Fig. m). 


Next we consider the momentum distribution which is 
an observable to gain information about cold atomic gas 
systems. The momentum distributions of the impurity 
for gBB = 0 and ujab = 1 is shown in Fig. |4]for differ¬ 
ent mass ratios and interaction strengths. These distri¬ 
butions can be understood from the discussions above. 
When the mass ratio increases the impurity wave func- 


As a final characteristic of the Bose polaron, we 
consider the overlap between the non-interacting and 
strongly interacting states for different values oi gBB and 
mass ratios as function of Nb- This quantity is related to 
the orthogonality catastrophe Q and has generated re¬ 
cent interest as a probe of many-body physics with cold 
atoms [H, [bll, In Fig. [5] we see a power-law behav¬ 
ior, but more interestingly, the exponent changes with 
both gsB and mass ratio. The overlaps remain finite for 
finite system sizes and only goes to zero for iV^ —>■ 00 
^|. The original work of Anderson Q uses a potential 
to model the impurity which corresponds to the limit 
mAB 00 . Consistent with Anderson, this limit shows 
very fast decay (high negative power dependence on Nb) 
but already for mass ratio mAB = 3 the suppression is 
considerable as seen in Fig. [5l On the contrary, in the 
opposite limit of equal masses we see much longer tails. 
Experiments using equal mass two-component setups and 
two atomic species with different masses could therefore 
complement each other perfectly when studying the or¬ 
thogonality catastrophe for Bose polarons. 


Experiments. Our predictions should be addressable 
using current experimental setups. In particular, effec¬ 
tive ID systems have been produced that exhibit be¬ 
havior consistent with zero-temperature predictions for 
both bosonic [l3 - [l^ and fermionic atoms Two- 

component bosonic systems in ID S [iinei can be 
used to explore the equal mass Bose polarons. Mass- 
imbalanced Bose-Bose mixtures in ID have been explored 
with ®^Rb and {mAB < 1) [1^ and new experiments 
with ®^Rb and ^^^Cs {mAB > 1) appear promising if an 
effective ID geometry can be reached 23 . Our theory 
provides predictions for experiments in the ID regime 
taking into account any experimental features such as 
different trap frequencies for different atoms, relative dis¬ 
placement of the trap, and mass imbalance. 
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FIG. 5: (color online) The overlap between the total wave 
function of the system aX g — Q and g = 50, | (V'g=o|'!/' 9 = 50>P, 
as a function of Nb- The best £t and the corresponding pa¬ 
rameters are shown in the figure. The upper value of {j3, a) 
corresponds to the upper curve, etc. For the infinite mass im¬ 
purity case, the fit parameters are qualitative as there is larger 
uncertainty here due to the small numerical values involved. 
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I 

Supplemental Material 


The model 


In the main text we consider an impurity particle interacting with Nb bosons. There we introduce a decomposition 
of the total wave function in the following form 

.. ,yNB) = '^(l}j{x)‘^j{yi, ■ ■ ■ ,yNB\x), (4) 

where x {yi) is the coordinate of the impurity (majority) atom, and ‘hj is a normalized eigenstate of the eigenproblem 
'Yl!i=i^B{yi)^j = Ej{x)^j where we assume that there is a zero-range delta-function interaction sitting at a fixed 
value of X. Thus we solve the eigenproblem as a function of the parameter x. Since there are no intra-species 
interactions among the B particles, we assume a product wave function for the bosons, i.e. 


Nb Nb 

■ ■ ^VNbU) = sWf^{yi\x), Ej{x) = Y^eyXx), (5) 

i=l * i=l 

where S denotes the symmetrization operator acting on the yi coordinates, and ff^i{yi\x) is the fc^th normalized 
eigenstate of the one-boson problem for fixed x with energy {x). Since the position of the impurity is fixed we use 
free solutions of the corresponding one-body problem on both sides of x and connect those with the delta-function 
boundary condition, i.e. 


f{x\x) = ^ 


df{y\x) 

dy 


y=x+E 


df{y\x) 

dy 


( 6 ) 


where the subscript y = x + e{y = x — e) denotes the right (left) derivative at point x. For convenience we drop 
index ki in and i in yi. Notice that when 1/g = 0 the boundary condition simplifies to f{x\x) = 0. As in the 
main text, we assume a harmonic confinement and adopt natural harmonic oscillator units for a majority particle. 
Free solutions for the harmonic oscillator potential without requirements of continuity and square-integrability are the 
Tricomi and Kummer confluent hypergeometric functions multiplied by a Gaussian function (7 (tI | . For our problem we 
need to choose such combinations of these free solutions such that the resulting wave functions are square-integrable 
and continuous together with their first derivatives everywhere except y = x. Moreover at y = x the wave functions 
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should satisfy Eq. For the sake of argument let us consider a; < 0, since the other case is obtained trivially after 
reflection. We obtain the following wave function for y < x 


f{y\x) = A{x)e 


(7) 


where a{x) = 1/4 — e{x)/2, A{x) is a normalization factor, U and M are the Tricomi and Kummer confluent hyper¬ 
geometric functions respectively. The solution for x < y < 0 reads 




f{y\x) = Aixje ■ {-V {/ 


And on the last interval, i.e. y > 0, we have 


f{y\x) = A{x)e • U ^a(a;), ^ 


-1 + 2 


M (a(x), i,0) 

U (a(x), 1,0) M (a(x), x^) 


M (a(x), O) C/ (a(x), |, x^) 
The boundary condition from Eq. m leads to the following equation for a(a;) 


( 8 ) 


(9) 


—2xa(x) U(a(x),l,0) / M (a{x), 1, x'^) 

g M(a(a:),i,0) {u{a{x),l,x‘^) 


U a{x) + l,-,x'^] +2M [ a{x) + 


+ 


U a{x), -,x'^] - 2 


U {a{x), ^,0) 
'M {a{x), i,0) 


M a{x), -,x'^]= 0. 


( 10 ) 


It is interesting to note that if the impurity is placed at the origin, i.e. x = 0, the above condition can be simplified 
using the asymptotic values: lim 2 ._,.Q+{7(a(x), i, x^) = Y(i/'^a{o)) l™a;->o+77'(a(x), i,x^) = ■ With 

these values Eq. OT yields the equation for the energy of two interacting atoms in a harmonic oscillator |65[ 


( 10^ ^ r(5/4-e(0)/2) _ 

(e(0) / )+^^r(3/2-e(0)/2) 

Once the functions f{y\x) and e(x) are determined we write down a system of equations for (j)j 


( 11 ) 


1 


rriAB 


i=i 


[Ha{x)+ Ei{x)]cl)i = , Pij{x) = Qij{x) = ( 12 ) 


dx 


d 


1 


52 


dx 


dx^ 


where itiab is the ratio of masses of the impurity and a boson, the subscript y on the brackets denotes integration 
over all yi,..., 2 /jvb- Since the function / is normalized for every value of x we derive that Pa = 0 and Qa < 0 571. 


Let us now estimate the relative contribution of the couplings terms. We start with Pjjt , this term is not zero if 
‘hj contains only one excitation of bosons compared to , i.e. if the corresponding products in Eq. ([S]) differ by 
only one function /. For bosons, excitation is a coherent process. We thus expect Pij to grow with Nb as \/Nb 
for the ground state. To understand this in detail we adopt the standard notation of second quantization, i.e. we 
write ‘hj = |ni, 712, n^), where n \,..., rik are occupation numbers in 1,..., k states. Without loss of generality let us 

assume that = \ni -I-1,712 — 1, •■•,7ife). The coupling is then 


/ d 

f2{y\x)^fi{y\x)dy, 


(13) 


where is the creation operator. This means that Pjj' ~ \/ 712(711 + 1), so we see that it is smallest if ni = 0, 772 = 1 
and largest if both tii ^ 712 ^ Nb/2. Similarly Qij is non-zero between states with either zero or one or two relative 
excitations. Obviously Qu ~ Nb, whereas if the excitation is only one then in analogy with Pij we have Qij ^ 
\/7i2(7ii + 1), where 711,712 have the same meaning as above. For two excitations let us consider only the following 
states (tii, 712, 7ife| and {ni -I- 1, 77,2 + 1, 773 — 1, 774 — 1,..., nfc| for which we obtain that Qij ~ -^ 773774(771 -I- 1)(772 -I- 1). 
Hence we conclude that only terms Qij that correspond to none or two excitations should be kept for large Nb A one 
considers the ground or low-lying excited states. From physical intuition we expect the coupling terms to contribute 
only if the impurity is in the middle of the trap, which corresponds to large values of At the same time, we expect 
that for strong interactions the impurity is pushed to the side of the trap to minimize the total energy of the system 
















in the ground state and in low-lying excited states. This allows us to neglect off-diagonal coupling terms Qy ■ Let us 
for clarity consider the ground state, then we find the impurity wave function, tpi (x), numerically from the equation 


Ha{x) + Nsegsix) 


Nb / f dfgs{y\x) 


TTlAB \ V 


Ui(a:) = E(j)i{x) 


(14) 


where fgs{y\x) is the ground state wave function of the one-boson problem for fixed x with energy tgs{x). Notice that 


the energy E provides a variational upper bound to the exact energy |56|. 


The density and pair-correlation function. 

In our analytical results we use the A^-body wave function, 'tj;{x,yi,... ^yNs) to obtain the densities of impurity 
and majority components, 


n{x) = J \iljfdyidy 2 ...dyNB, nB{y) = Nb j |'0|^dxd?/2dy3 ■ • ■ dyjvB, 

and the pair-correlation function for a boson-impurity pair, 

nAB{x,yi) = j \ip\'^dy2 ■ ■ - dyNg. 

To find the momentum distribution we first define the Fourier transform of the wavefunction: 

1 -\-Nb p 


(15) 


(16) 


V'(p,gi,.. ■,qNB) = 


1 




/ all space 


V'(xi,j/1,.. .,yNBy^^e^‘’^y^ . ..Ey^BV^Bdxidyi. ..dyNB, (17) 


where p {qi) is the momentum of the impurity (majority) particle. Next we obtain the momentum distribution for 
the impurity using the assumption that the wave function is factorized. 


{p) = [ (t>*{x)(j){x) ( [ f*{y\x)f{y\x)] dxdx. 

J all space \J / 


(18) 


The Gross-Pitaevskii equation in a quasi-one dimensional geometry 

For convenience of the reader we review the Gross-Pitaevskii equation (GPE) which describes the ground state for a 
system of weakly-interacting identical bosons at zero temperature that is specified with the following three-dimensional 
Hamiltonian with h = 1 


Nb 




2m B 


'■3)^ 


(19) 


i<i 


where is the coordinate of the Tth boson, tob is the mass of a boson. For small values of the boson-boson scattering 
length as the coupling constant reads Gbb = see e.g. [Tlj. The external potential is taken to be a harmonic 

oscillator, V{ri) = ^m{uj‘^x^ +(jJyy^ -I-w^z^). Finally, i5(r) is the 3D Dirac delta-function. This Hamiltonian describes 
a weakly-interacting system meaning that the s-wave scattering length is much smaller than the average distance 
between atoms. To obtain the GPE we assume that the total wave-function 41 of the system of Nb bosons is a 
product of single-particle functions tp, 


4'(ri,r2,..., vnb) = ipiri)'ipir2). ..ipivNB), 


( 20 ) 


which minimizes the expectation value of the Hamiltonian. To satisfy this requirement the single-particle wave- 
function should solve the GPE, 


+ '^^( d ) + GbbNb\P^\^ ) 'ip, 


(21) 
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where /i is the chemical potential. The chemical potential appears in the equation as a Lagrange multiplier from the 
normalization condition for the wave function, i.e. f dV\tl’\'^ = 1. To simulate one-dimensional geometry experimen¬ 
tally one confines a system by a potential with two frequencies that are much larger than the third. At very low 
temperature this means that any two atoms outside of the interaction range move effectively in ID. Having this in 
mind we scale the GPE using the units given by the frequency of weak confinement: r —>■ r6 and t) —>• 

where h=^ —-—, such that 

Y rriBi^y ’ 






+ H(r) 


AttUsNe 






( 22 ) 


where = LOxju^y and 72 = Let us now assume a cigar-shaped trap with ^ Wy. This means that 

the time evolution does not cause excitations along the x- and 2 ;-axis, since these require much more energy than 
the excitations along the y-axis, see e.g. [s^. We therefore assume that a single-particle solution of the GPE has 
the following form: y, z) = ' 4 ’\{y)ip 2 'i{x, z), where the function in x- and z-coordinates is the ground state of the 

harmonic oscillator: ■023 = This means that {x,z) part can be easily integrated out which 

leads to 


where gss is given as: 


^ / 

—fpi = 

UJy 


1 


1 


“2 ^ ”^2^ + A^b ■5sb|^/’i| ]^i 


AtTUs 
9bb = -r- 


\il^2‘i\^dxdz 


A-Kas y/lTh 
h 2 tt 


^Qjsy^uJxdJz 

buiy 


(23) 


This equation supplemented with the delta-function boundary condition at the point where a boson meets the impurity 
is used in the main text to describe weak interactions among the majority particles. 










